Lowest eigenvalues of the Dirac operator for two color QCD at nonzero 
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We investigate the eigenvalue spectrum of the staggered Dirac matrix in SU(3) and U(l) gauge theory as 
well as in full QCD with two colors and finite chemical potential. Along the strong-coupling axis up to the phase 
transition, the low-lying Dirac spectrum of these quantum field theories is well described by random matrix theory 
and exhibits universal behavior. Related results for gauge theories with minimal coupling are discussed in the 
chirally symmetric phase and no universality is seen for the microscopic spectral densities. 



1. Spectrum in confinement 

The eigenvalues of the Dirac operator are of 
great interest for the universality of important 
features of QCD and QED. The accumulation of 
small eigenvalues is, via the Banks-Casher for- 
mula [0, related to the spontaneous breaking of 
chiral symmetry. Their properties are known to 
be described by random matrix theory (RMT) in 
the confinement, see Ref. 

The Dirac operator itself, p = ^ + ig^, is anti- 
Hermitian so that the eigenvalues A„ of ip are 
real. Because of {.^,75} = the nonzero A„ oc- 
cur in pairs of opposite sign. In the presence of 
a chemical potential > 0, the Euclidean action 
becomes a complex number and the fermionic ma- 
trix becomes a non-Hermitian matrix. 

2. Spectrum into deconfinement 

We have continued our investigations |^ with 
a study of the distribution of the small eigen- 
values in the whole phase diagram. The Banks- 
Casher formula ^ relates the Dirac eigenvalue 
density p(A) at A = to the chiral conden- 
sate, E = \{i>il>)\ = linie^o limy^cx) 7rp(e)/V. 
The microscopic spectral density, Ps{z) = 
limv^oo p{z/VT,) /VE, should be given by the 
appropriate prediction of RMT Q . In the case of 
complex eigenvalues the situation is more compli- 
cated 1^ and p(0) can be used as a lower bound 



for E. 

We present results for the density of the small 
eigenvalues in Fig. Q for SU(3) theory and the 
staggered Dirac operator on a 4* lattice from 5000 
configurations for (3 = 5.4 and 3000 configura- 
tions for (3 — 5.6 and f3 — 5.8 around the crit- 
ical temperature Pc ~ 5-7- Our analogous pre- 
sentation for U(l) theory is from 10000 config- 
urations on a d'* lattice around the critical cou- 
pling (3c ~ 1.01. In the confinement phase it was 
shown that, both the microscopic spectral density 
Ps{z) and the distribution P(Ainin) of the small- 
est eigenvalue agree with the RMT predictions of 
the chiral unitary ensemble for topological charge 
ly ^ [Hj^. In the case of two-color QCD with 
staggered fermions on a G'* lattice we produced 
at least 2100 configurations for each value of /i 
around the critical chemical potential pc ~ 0.3 
keeping ^ = 1.3 fixed §. Since the eigenval- 
ues move into the complex plane for /i > 0, a 
band of width e = 0.015 parallel to the imagi- 
nary axis was considered to construct p{y), i.e. 
p{y) = dx p{x, y), where p{x, y) is the density 
of the complex eigenvalues x + iy. Alternatively 
p(|A|) was constructed from the absolute value for 
|A| small and is presented in this write-up. 

The distribution of the lowest eigenvalue is dis- 
played in Fig. H for the theories under investiga- 
tion. The data of the SU(3) and U(l) theories 
in the deconfinement could not be fitted to the 
width of the functional form derived for RMT in 
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Figure 1. Density p{\) of small eigenvalues for SU(3) (left) and U(l) gauge theory (center) as well 
as for two-color QCD (right) across the transition of critical coupling and critical chemical potential, 
respectively. 
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Figure 2. Distribution P(A,nin) of the smallest eigenvalue as in Fig. |^. 
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Figure 3. Fit of the spectral density to p(A) = C(A - Ao)2™+i/2 in SU(3) at /3 = 5.8 (left), in U(l) at 
P = 1.10 (center) and in two-color QCD at /i = 0.4 (right). The contribution of the smallest eigenvalue, 
the 11*^ eigenvalue and the 21^* eigenvalue is inserted. 
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the chirally broken phase. The quahty of the data 
of two-color QCD was not sufficient for a rehable 
fit to a trial function for P(|Ainin|) except for the 
RMT prediction of the chiral symplectic ensemble 
at /i = 0. The lattice results of the three theories 
look very similar and might give some help for a 
derivation within RMT. 

Nevertheless, the quasi-zero modes which are 
responsible for the chiral condensate S ^ build 
up when we cross from the deconfinement into the 
confined phase. Figures |l| and |^ demonstrate that 
both p{X) and P(Aniin) plotted with varying /? or 
jj, on identical scales, respectively, can serve as an 
indicator for the phase transition. The change 
in bending from positive to negative curvature of 
p(0) as a function of coupling/density might serve 
to pin down the critical point. 

In Fig. ^ we turn to a discussion of the spec- 
trum in the quark-gluon plasma and Coulomb 
phase. From RMT a functional form of p(A) = 
C(A - Ao)^™"*""^/^ is expected at the onset of the 
eigenvalue density [||. A fit to the data in the 
regime up to A = 1 yields m — 0.081(1) for 
SU(3) and m = 0.084(1) for the Abehan the- 
ory, in agreement with recent studies This 
suggests that both theories correspond to univer- 
sality class TO = 0. For this class a microscopic 
level density involving the Airy function can be 
deduced from RMT |0 . A rescaling of our data 
from the 4'^ lattice to this functional form is not 
satisfactory for both theories |9). We tried a fit 
for two-color QCD assuming the above functional 
form for p(|A|) in a region up to |A| = 1 obtain- 
ing TO = 0.03(1) consistent with universality class 
TO = 0. 

3. Conclusion 

We investigated universality concerning the 
low-lying spectra of the Dirac operators of both 
SU(3) and U(l) gauge theories. One finds that 
in the phase in which chiral symmetry is sponta- 
neously broken the distribution P(Amin) and the 
microscopic spectral density Ps{z) are described 
by chiral RMT. For two-color QCD with chem- 
ical potential we were not able to verify such a 
relation for < /i < /^c, not only because of our 
data but also in lack of an analytic result for non- 



Hermitian RMT ||ll|. In the phase where chiral 
symmetry is restored one has to rely on ordinary 
RMT. Here we find universal behavior only of the 
macroscopic density p(A) for all gauge theories 
with minimal coupling. 
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